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Abstract 

This is the continuation of Montes' paper "On the canonical discussion of polynomial 
systems with parameters" . In this paper we define the Minimal Canonical Comprehen- 
sive Grobner System (MCCGS) of a parametric ideal and fix under which hypothesis 
it exists and is computable. An algorithm to obtain a canonical description of the 
segments of the MCCGS is given, completing so the whole MCCGS algorithm (imple- 
mented in Maple). We show its high utility for applications, like automatic theorem 
proving and discovering, and compare it with other existing methods. A way to detect 
a counterexample is outlined, although the high number of tests done give evidence of 
the existence of the MCCGS. 

Keywords: comprehensive Grobner system, canonical, minimal, reduced specifica- 
tion, generalized canonical specification, constructible sets. 
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1 Introduction 



In this paper we continue the task introduced in dMontes 200iD . Let us briefly remember 
the basic features. 

Given a parametric polynomial ideal / C -fir[a][x] in the variables x = {xi, . . . , Xn) and 
the parameters a = (ai, . . . , am), and monomial order y^, our interest is to find the different 
types of solutions for the different values of the parameters. Let i^' be a computable field and 
K an algebraically closed extension. A specialization is the homomorphism : i^[a][x] — > 
that corresponds to the substitution of the parameters by concrete values a £ K . 
A comprehensive Grobner system (CGS) is a set of pairs: 

CGS(/, y^) = {{Si,Bi) : Si C ;r^constructible sets, Bi C A[x], 

aa{Bi) = gb((Ta(/), ^^r) Va G 5i , and jj. Si = 1^}, 



*This research was partly supported by the Ministerio de Ciencia y Tecnologia under project MTM 
2006-01267, and by the Generalitat de Catalunya under project 2005 SGR 00692. 
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where the Si are caUed "segments" and the Bi "bases" . Frequently the word "segment" is 
also used for the pair {Bi, Si) whenever the sense is clear from the context. 

There are different known algorithms that provide Comprehensive Grobner Bases and 
Systems for a given ideal: 

CGB (jWeispfenning 1992l l. 

ACGB dSato-Suzuki 20031 : ISato 20051 : iNabeshima 2005l 'l. 

SACGB dSuzuki-Sato 20061 ^. 

HSGB (|Gonzalez-Traverso-Zanoni 20051 ). 

BUILDTREE (|Montes 20021 : iManubens-Montes 20061 : iMontes 200l\ ). 

There are available implementations (jPolzman-Seidl-Sturm 20061 ) of W eispfenning's CGB 
algorithm in Reduce, of Suzuki-Sato's SACGB in Risa/Asir and in Maple ( Suzuki-Sato 20061 ) 
and of Montes's BUILDTREE in MapleS. All these algorithms allow to build both Com- 
prehensive Grobner Bases and Systems, but they are differently oriented. A comparison of 
the most interesting among them is given in section [6l 

In fact, comprehensive Grobner systems are in general more effective to handle for their 
use in the applications than comprehensive Grobner bases. But it is also convenient to 
require some more additional features to these Grobner systems when looking for applica- 
tions. 

The first requirement is to have disjoint and reduced CGS. By disjoint we mean that 
the Si form a partition of K , and by reduced that the bases Bi specialize to the reduced 
Grobner basis of <Ta(/) preserving the leading power products (Ipp), for every value a 
of the parameters inside Si. The algorithm BUILDT REE (introduced in ( Montes 20021 1 
as DISPGB and improved in ( Manubens-Montes 20061 )) already builds a disjoint, reduced 
CGS. 

In (jMontes 2007l l the interest is focused in the improvement of BUILDTREE to obtain 
a simpler and canonical CGS. The method consists of grouping together all the segments 
with the same Ipp that allow a same basis specializing well on all the grouped segments. A 
natural conjecture establishes the existence of an equivalence relation between the segments 
having the same Ipp, and an algorithm is given to compute the basis corresponding to the 
grouped segments. 

In order to obtain a truly canonical CGS we nee d to describe t he segments in a canonical 
way. This is the objective of the present paper. In (jMontes 20071 ) a canonical description of 
a segment determined by a diff-specification was already given, but it remained to obtain a 
canonical representation of the addition of such segments. The objective is thus to obtain 
the MCCGS (minimal canonical CGS). 

Definition 1. We cah Minimal Canonical CGS (MCCGS) a CGS with the following prop- 
erties: 

i) disjoint CGS, i.e. Si n = for i ^ j; 



^The library DPGB 7.0 written in Maple 8 is available at the web http://www-ma2.upc. edu/'^montes 
and is actualized with the MCCGS algorithm. 
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ii) reduced CGS, i.e. the polynomials in Bi have content 1 w.r.t. x, Bi specializes to the 
reduced Grobner basis of (Jail) for every a G S'j, their leading coefficients are non-null 
on Si and their Ipp remain stable; 

iii) the sets Si, . . . ,Ss are intrinsic for the given / and >-x and are described in a canonical 
form. 

iv) the number of segments of the CGS with the above properties is minimal. 

The currently existing algorithms that can build comprehensive Grobner systems, say 
BUILDTREE, CGB, CCGB, ACGB and SACGB, do not hold all these properties. BUILDTREE 
builds a comprehensive Grobner system satisfying properties i) and ii). But CGB, ACGB 
and SACGB do not hold property i), at least. Finally, although the Grobner system ob- 
tained within CCGB is canonically determined, does not hold properties i) nor ii) as for 
the obtention of a comprehensive Grobner basis the algorithm needs the Grobner systems 
to be faithful. 

It must be emphasized that the existence of the MCCGS depends on the Conjecture for- 



mulated in (jMontes 20071 ) about the existence of an equivalence relation between segments 
allowing a common basis. 

If the Conje cture is true, then the computation using MCCGS algorithm proposed 



in (jMontes 20071 ) and in this paper, already depends on the semi-algorithm GENIMAGE 
given there for computing pre-images, that uses arbitrary bounds. 

With these restrictions, MCCGS algorithm builds a comprehensive Grobner system 
satisfying all the properties in Definition [TJ These properties will make the algorithm more 
suitable for the applications. In particular, they are very appropriate for automatic theorem 
proving and discovery (see ( Montes-Recio 200?! )) as well as to compute geometric loci as 
shown in example [H 

Furthermore, MCCGS also allows to restrict the parameter space to a constructible 
set and impose a-priori null and non-null conditions. This is also an interesting tool for 
applications when we want some degenerate cases to be avoided (see Section [5]) or restrictions 
on the parameters to be given. For example, when the parameters involve angles, and the 
equations are given using the sine and cosine of the angles as parameters, it is important 
to restrict the solutions to cos^ 99 + sin^ (/j — 1 = 0. 

The whole algorithm MCCGS is achieved by three steps: 



BUILDTREE (described in (jManubens-Montes 200fil )) 



grouping segments with common basis (described in ( Montes 20071 )) 



iii) representing the subsets in canonical form. This part will be described in sections [3] 
andU 

Although the algorithm requires two term orders (one for the variables and another for 
the parameters >-a)-, the result will not depend on >-a, as the segments {Bi, Si) are intrinsic 
for the given ideal / and the term order >-x- Even though, >—a will be used to determine 
the reduced Grobner bases of the ideals involved in the description of Si. 

The pape r is structured as follows: section [2] is devoted to recalling some properties and 
results from ( Montes 200?! ) which are used in the subsequent sections. The generalization 



3 



of the canonical specification and its properties are given in Section [3l In Section [J] we 
give the algorithm which collects the corresponding segments into a generalized canonical 
specification and builds up the Minimal Canonical Comprehensive Grobner System (MC- 
CGS). In section [5] a practical application to automatic theorem proving is given. Finally, 
in section Owe compare the main available CGS algorithms. 



2 Preliminaries 

We describe now briefly steps i) and ii) of the MCCGS before tackling the last step iii) that 
will be studied in this paper. The algorithm starts with a parametric ideal / and a term- 
order >—x on the variables. An auxiliary term-order >—a over the parameters is needed to 
describe the subsets in K using Grobner bases. It does not affect the segments themselves 
but only their description. 

Step i) is per formed by BUILD TREE algorithm, and w as described for the first time 
in ( Montes 2002 ) and improved in ( Manubens-Montes 20061 ) . The output is a disjoint re- 



duced CGS, where the subsets Si are determined by red- specifications. A red-specification 
of a segment S is described by the pair {N, W), where N is the radical null-conditions ideal, 
and is a set of irreducible (prime) polynomials on K\a] representing non-null conditions 
such that no prime component Ni of the prime decomposition of N does contain any of the 
polynomials in W . We have S = N{N)\Y{h) with h = ritoGH/ ^- ^ red-specification deter- 
mined by (A^, W) is easily transformed into a diff- specification (N, M) with A^ C M where 
S = V(A) \ V(M), by considering the polynomial h = Hwety ^ ^^^^ taking M = (h) + N . 

Let us denote CGSi the output of BUILDTREE that consists of a list of segments each 
represented by the three objects {Bi, Ni, Wi). Remember that each of these segments have 
characteristic set of Ipp of their bases Bi that are preserved by specialization on Si. We 
say that a basis G specializes well to {B,N,W), with Ipp(G) = lpp(i?), if the polynomials 

— N 

of G are proportional to the polynomials of B, i.e. for each g £ G there exist f £ B and 
a,l3 eW* such that ag'^ = /3 f, where W* = {k l\t=i wf' ■.keK,Xie Z>o, Wi G W}. 



Step ii), described in ( Montes 200?! ). selects the segments of CGSi with the same Ipp 



that admit a common reduced basis specializing well to the r educed Grobn er basis for 
every specialization in the grouped segments. If Conjecture 7 in (jMontes 2007l l is true, the 



grouped segments form an intrinsic partition of the parameter space. To perform that task 
the algorithms DECIDE and GENIMAGE are used. The first one tests whether one from 
two segments with the same Ipp has already a generic basis specializing to the other (this 
is the most frequent case) or a sheaf exists and is necessary or whether possibly a more 
generic basis must be found (by GENIMAGE). Whenever no pre-image nor sheaf is found 
then both segments are not equivalent and cannot be summarized. It can happen that 
instead of simple polynomials the basis Bi contains also sheaves of polynomials. A sheaf 
{ffi) • • • ,9k} is accepted in a basis of a segment instead of a simple polynomial, whenever 
all the polynomials in the sheaf specialize to the corresponding polynomial of the reduced 
Grobner basis of the specialized ideal or to 0, and some of the polyno mials in the sheaf 
specialize to non-zero for every Tx £ Si. As was shown in (jWibmer 2006l l. it is necessary to 



use sheaves for some over-determined systems if we want to group all the segments admitting 
a common basis with the same 1pp. We must notice that DECIDE algorithm also depends 
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on the semi- algorithm GENIMAGE to determine a polynomial / that specializes well to fi 
over (A''i, Wi) and to /2 over (A^2, ^^2)- Thus the canonicity of the results of the computation 
of a MCCGS relies on GENIMAGE and the truthfulness of the mentioned conjecture. 

Let us denote the output of the second step CGS2. It will be described by segments 
with a common basis i?j and a set of red-specifications: 

{Bi,{{Na,Wa),...,m,^,W,,J}). (1) 

Si will now be the union of the segments determined by the red-specifications (A^ifc, Wj/c) 
for k from 1 to ji . 

Step iii) will be described in next sections. Its ob jective is to gi ve a canonical description 



of the union of the grouped segments of step ii). In (jMontes 20071 ) it was shown how a diff- 



specification can be transformed into a can- specification. Here we will prove that the union 
of red-specifications or their corresponding diff-specifications can be transformed into a 
ge neralized can-s pecification using what we call a P-tree. The idea is based on Theorem 12 



in (|Montes 20071 ). Let us give here a slightly different formulation of it, more appropriate 



for the current purposes. 
Theorem 2. 

i) Every diff- specification S = Y{N) \ V(M) admits a unique can- specification 

S = Y{N) \ V(M) = U (V(iV,) \ (U,V(M,,-)) , (2) 

i 

where M = DiNi and Mi = DjMij are the irredundant prime decompositions over A 
of the radical ideals M and Mi respectively, where Ni C Mij. 

ii) The Zariski closure over IV^ verifies 



S = [j {V{Ni) \ (U,-V(M,,)) = U V(iVi) = V(AA). 

i i 

iii) The can- specification verifies 

V(iV,) \ (UjV(Mij)) = S^^N{Ni). 



iv) Given a diff- specification of S the algorithm DIFFTOCANSPEC (jMontes 20071 ) builds 
its can-specification. 

The need of having a canonical description of the intrinsic segments comes from the 
need of comparing different outputs for the same problem, and also to have a final simple 
description of the segments. 
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3 Adding segments 



We tackle now the third step of MCCGS, i.e. the description of the union of the segments in 
a canonical form. We start with segments of the form ([1]). The red-specifications {N, W) can 
be transformed into diff-specifications {N,A4), as explained in Section [21 so we are attained 
with the obtention of a canonical representation for the addition of diff-specifications. We 
cannot assume that the simple form given by formula ([2]) will be sufficient. A more complex 

constructible set will be formed grouping all the segments Sik for 1 < k < .?,;. 

Thus we generalize the concept of canonical specification given in ([Montes 20071 ): 

Definition 3 (P-tree). A P-tree is a rooted directed tree such that 

i) the nodes are prime ideals over A except the root, denoted r, 

ii) when P ^ Q is an arc then P Q, 

iii) the children of a node are a set of irredundant prime ideals over A, (whose intersection 
form a radical ideal). 

By definition the root level is 0. 

Definition 4 (C-tree). To any P-tree we associate an isomorphic C-tree by changing every 
node P to a subset of K denoted C{P) by the following recursive procedure: 

i) if P is a leaf (terminal vertex) then C(P) = V(P), 

ii) if P is an inner node different from the root and Pi , . . . , P^ are its children, then 

C{P)=Y{P)\{C{Pi)U---UC{Pd)) (3) 



iii) if Pi , . . . , P(i are the children of the root vertex r then 

C{r) = C{Pi)U---UC{Pd). 



Note that for C(r) the parity of the vertex-level acts additively for odd level vertices 
and as a subtraction for even level vertices. (See example [6] below) . 

Definition 5 (Generalized canonical specification). A generalized canonical specification 
(GCS) of a set S* is a P-tree such that S = C{r) satisfying, for every node P at level j, the 
following condition: 

C{P) = V(P) n B (4) 
where P = for j odd and B = A™ \ S for j even. 

Example 6. To clarify the definition suppose that we want to describe the set 5*1 of the 
M^-space with coordinates a, b, c consisting of the planes a = and b = —1 except the lines 
a = b = and a = c = plus the point 0(0, 0, 0). We can express Si as 

5i = ((V(a) U V(6 + 1)) \ (V(a, b) U V(a, c))) U V(a, b, c) 
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Si 52 




[a,b,c] [a,b,c] [a,b + l,c] [a,b,c] [a,b,c] 



Figure 1: Trees representing the sets and 52 in generalized can-specification. 



But there exist many other possible determinations of this set. If we want to obtain the 
GCS of Si we must write Si in the form 

Si = (V(a) \ ((V(a, b) \ V(a, b, c)) U (V(a, c) \ (V(a, b, c) U V(a, b + 1, c))))) U V(6 + 1) 

This formula can be represented by the tree associated to 5*1 shown in figure [TJ Notice that 
we must include the point V(a, 6+ 1, c) under the branch of ¥(a, c), as this point belongs to 
5i and condition @ requires it to belong also to C (a) . The interest of that representation 
lies in the fact that it is unique as we prove in Theorem [7] below. 

Consider now the set S2 = Si\ V(a, 6+1, c). In order to preserve property @ of the 
GCS definition, the P-tree associated to S2 will be modified from the P-tree associated to Si 
by eliminating the point under the variety V(a, c) and setting it under the variety V{b + 1). 
The new tree is also shown in Figure [TJ These examples should clarify the definition of 
GCS to obtain canonicity of the description by preserving condition 

Theorem 7. A subset S C defined by a GCS has the following properties: 

i) For every vertex P, except for the root, 

C{P) = V(P) 

where, as usual, the Zariski closure is taken over 'K^ . 

ii) For the root vertex r 

d 

1=1 

where the Pi 's are the children vertices of r. 

iii) S has a unique GCS decomposition. 

Proof. i) The inclusion C is obvious as C(P) C V(P). To prove the equality we have 

d d 
C{P) = V(P) \ IJ C{P{) D V(P) \ U V{Pi). 

1 = 1 4 = 1 
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Consider the closure of the above formula and apply Theorem [2] (ii). The result 
follows. 

ii) Is an immediate consequence of i). 

iii) To prove the uniqueness we proceed by induction on d. For d = 1, the tree is formed 
by the root r and a set of children nodes forming an irredundant prime decomposition 
of the radical ideal defining S, by Definition [3] iii) . Thus, in this segment the P-tree 
is unique. 

Assume now by induction hypothesis the uniqueness of the GCS for every P-tree of 
maximum depth less than d and let us prove, as a consequence, the uniqueness also 
for depth d. Let S be defined by a P-tree of maximal depth d representing a GCS. 
By part (ii) of the Theorem we have 

i 

where the Pj's form the unique irredundant prime decomposition over A of the rad- 
ical ideal rijPj defining S by Definition [3] (iii). Thus they are uniquely determined. 
Denoting Pij the children of P,, by (jH), we have 

di 

C(P,) = V(P,) \ U C{P,,) = Y{Pi) n S (5) 
i=i 

showing that C{Pi) is also uniquely determined. Set Si for the subtracting set 

Si = [JC{P,,). (6) 
i=i 

As Si C V(Pi), Si is also uniquely defined by By Definition [5l formula (0]), we 
have 

CiPij)=Y{Pij)n{K"'\S), 

Thus 

s^ = [J C{P.,) = I [J v(p,,) I n (Z™ \ s) 

and so 

CiPi,)=YiPi,)nSi (7) 

By the ascending chain condition for the ideals in the branches and condition ^ for 
the P-tree of S, equation ([7|) ensures that condition ^ is also respected for the subtree 
of Si, whose root vertex is given by ([6]). Thus the subtree of Si also forms a GCS of 
Si with depth less than d. By the induction hypothesis it is uniquely determined and 
so does the complete P-tree of S. 

□ 
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4 The MCCGS algorithm 



Given an ideal / and the monomial orders >-x for the variables and >—a for the parameters, 
the following sequence of algorithms build up the P-tree T corresponding to the Minimal 
Canonical Comprehensive Grobner System associated to / and >-x- We describe them in 
descendent design. 

tree T ^ MCCGS(5, ^a) 

Input: B a basis of the parametric polynomial ideal I and monomial orders >-x, >—a- 
Output: T a tree containing the minimal canonical comprehensive Grobner system 
associated to I. 



To :=BUILDTREE(5, >-x, ^a) 
S :=SELECTCASES(ro) 
T :=GENCANTREE(S') 



MCCGS uses BUILDTREE (see (jMontes 20021 : iManubens-Montes 20061 )^ to build up 



the discussion tree Tq containing a CCS whose segments are expressed as red-specifi- 
cations. Then SELECTCASES takes Tq as input and classifies the segments from the 
CGS associated to Tq into pairs of the form where li is a set of red-specifications 

{{Nil, Wii), . . . {Nij^, Wij^)} whose corresponding bases have been generalized by the same 
basis Bi. Afterwards, MCCGS calls the new algorithm GENCANTREE to finally obtain 
the MCCGS associated to the initial ideal and term order. 

GENCANTREE uses GCS algorithm to build the P-tree corresponding to the general- 
ized canonical specification of the addition of segments. GCS algorithm begins by setting 
the ideal {0} at the root of new tree T and calls iteratively the recursive algorithm ADD- 
CASE. It must be noted that there are two kinds of nodes, namely odd level vertices and 
even level vertices, that are treated differently by ADDCASE. ADDCASE uses two aux- 
ihary algorithms: DIFFTOCANTREE (a minor transformation of DIFFTOCANSPEC) 
converts a diff-specification into a P-tree containing the associated can-specification, and 
SIMPLIFYSONS just makes the suitable simplifications. 

At the first iteration ADDCASE stores under root the P-tree of the unique canonical 
specification associated to {Nn, Wn). Then, to add each further red-specification {Nik, Wik), 
ADDCASE executes itself recurrently in post-order at the even level vertices u £ T and 
adds the can-specification associated to {Nik, Wik) contained in Y{Pu)- For example, in 
figure [2] it would act successively on the vertices 

c, /, i, j, A m, g, o, p, d, t, u, r, v, a. 



Thus, before acting on an even vertex u €z T, the algorithm must have acted on all its 
even descendants. Therefore, if an even level descendant v verifies that Nik i5 Pv, then the 
can-specification associated to {Nik, W^jfc) niust have been completely hung under v. In this 
case the test variable will contain false and thus DIFFTOCANTREE for current {Nik, Wik) 
will not act on P„ nor on any of its ascendant vertices. We must also remind that the 
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set of pairs S ^ SELECTCASES(ro) 

Input: To a BUILDTREE discussion tree whose terminal vertices shape a CGS with 

red-specifications. 

Output: S a finite set of pairs of the form {Bi, {{Nii,Wii), . . . , {Nij^,Wij.)}) taken from 
the CGS associated to Tq. 

G := {{Bi,Ni,Wi), iBr,Nr, Wr)} {the CGS associated to Tq} 

while G 7^ do 

Let {B, N, W) be the first element of G 
Bo := B; Nq := N; Wq := W; 
I := {{No,Wo)} 
G:=G\{{Bo,No,Wo)} 

for all {B',N', W) G G such that Ipp(-B) = lpp(5') do 

p:=0 

for all f ^ B while p 7^ false do 

Let /' G B' be such that lpp(/) = lpp(/') 

p := DECIDE(/, TV, W, f, N', W) 

if p 7^ false then 

Substitute / by p in Bq 

end if 
end for 

if p 7^ false then 

I :=/U{(iV',VF')} 

B := Bo; N ■.= NnN'; W := W n W; 
end if 
end for 

S := SU{{B,l)} 

G:=G\ {{B', N', W) G G such tliat (iV', W) G 1} 
end while 



ideals associated to the paths in T starting from root form ascending chains of prime ideals. 
Thus, whenever test is false, the condition cited above will also hold for all vertices placed 
between u and v, even the odd level ones, i.e. for all w; G T descendent of u and ascendant 
of V, Nik 2 Pw _ 

This way, ADDCASE completes current P-tree T to a new tree such that for every 
odd level vertex u with prime ideal P„, all points in V(P„) fl (V(iVjfe) \ V(/ijfe)) (where 
hik = I\weW,k ^) ™ C{Pu), as required. 

Nevertheless in the new tree completed by ADDCASE it could happen that P„ + N^k = 
Pu for some even level vertex u, which would cause that P„ and its unique child -Pchiid(w) 
coincide. If so, SIMPLIFYSONS takes the subtree under child(M), slips it upwards hanging 
it from parent(?/) and eliminates both vertices u and child(?i) from the tree. When this 
action is performed, it could also happen that some set of current even level siblings do not 
preserve the prime decomposition irredundancy property, as some lifted primes can contain 
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tree T ^ GENCANTREE(5) 

Input: S a finite set of pairs of the form {Bi, {{Nn, Wn), . . . , (iVjj., Wjj.)}). 
Output: the canonical tree T associated to S. 

initialize T 

for 1 < z < tlS do 

Create Ui a new vertex in T hanging from root 
store Bi in Ui 

I := {{Nii,Wii), . . . , {Niji,Wij^)} {red-specifications associated to Bi} 
T :=GCS(/) 
hang T from Uj 
end for 



tree T ^ GCS(Z) 

Input: I a finite set of red-specifications 

Output: a tree containing the Generalized Can-Specification associated to the addition of 
segments in I. 

initialize tree T with root r 

set Pr := 4> 

for all pairs (iV, W) G / do 

T :=ADDCASE((iV, W),r,T) 
end for 



some of their sibling vertices, i.e. 3vi,V2 G children(tt) such that Py^ C Py^ for u an even 
level vertex in T. SIMPLIFYSONS algorithm also detects these cases and eliminates the 
subtrees hanging from V2 as well as V2- Though, the action of SIMPLIFYSONS will restore 
the GCS-condition property of the tree. 

Note: For algorithmic reasons, all paths starting from root vertex in a P-tree will be of 
even length. Thus for odd length branches, the algorithm will add a new vertex [1] at the 
end. 

The above described algorithms build the complete MCCGS of the initial ideal. The 
following theorem states that GCS algorithm builds the generalized can-specification (GCS) 
associated to the set of the corresponding diff-specifications: 

Theorem 8. Given a finite list of pairs I = {{Nik, ^ife) ■ k = 1, . . . , M} of red- specifications, 
GCS(/) computes the P-tree associated to the generalized can-specification determining the 
constructive set 

M 

\JV{Nik)\Y{ n w). 

k=l w£Wik 

Proof Let S = Ufeli V(iVifc) \ ^iUweW.k'^)- The proof is done by induction on M, the 
number of red-specifications to be added. 

For M = 1, GCS uses DIFFTOCANTREE just once and, by Theorem 1 (iv), it builds 
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up the unique can-specification in tree T. Thus T is a P-tree such that C(T) = V(iVji) \ 

By induction hypothesis, assume now that after the M — I iteration of ADDCASE the 
GCS tree of the M — 1 red-specifications has been built and let T be this tree, which is a 
P-tree such that C(f ) = IJ^l^ ^{Nik) \ ^{UweWik ^) ^^^^ ^^^''^ vertex uef 
holds that C(P„) = V(P„) n C{f). We shall prove that the M-th iteration will build the 
GCS tree of S. 

Let us describe how the recursive ADDCASE algorithm acts on T adding Y(NiM) \ 
'^(YlweW M Denote by A(u) the operation on an even level vertex u that hangs to it the 
tree associated to the can-specification of {NiM, W^m) contained in V(Pu) (i.e. Y{NiM+Pu)\ 
V(A^iM + P?t + (n«)eVKM whenever it can be hung and returns false or true depending 
on whether parent(P) C NiM or not, respectively. So it hangs the points V(P„) fl (V(A^jM) \ 

A{u) is applied recursively in post-order. If A{u) returns false at some even level vertex 
u, the whole set Y{NiM) \ "^{NiM + (IlujeVK M ^^-^ been hung under u and thus, as u is 
even, C(father(M)) D V(A^jM) \ "^{^iM + (HtyGWiM ^))- Then A will not be applied to any 
of its ascendant vertices because C(T) = C{T) U ^V(A''jM) \ V(Arj^ -|- (HiueVKiM '"^))) • 

If A{u) returns true for all u which means that V(A'iM) \V(A^iM + (riiueH^M 
not completely been hung under root, then the P-tree corresponding to the red-specification 
{NiM, WiM) computed by DIFFTOCANTREE will be hung from root. Thus, we finally have 

that C(T) = C{f) U (YiNiM) \ Y{NiM + (R^ev^.M ^))) • 

This way, GCS algorithm obtains, as SIMPLIFYSONS ensures, a P-tree T such that 
for every node v eT holds that C(P^) = V(P^) n C(T) and C{T) =S. □ 
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(bool, tree T) ^ ADDCASE((iV, W),u, T) _ 
Input: (N, W) a red-specification , u the current vertex in P-tree T. 

Output: false if (A^, W) is not to be added to parent vertices, true otherwise. It also 
returns current tree T. 

if u is not terminal then 
test := treu 

for all V G children(u) do 
for all w G children(t>) do 

if AI)'DCASB{{N,W),w) = false then 

test := false 
end if 
end for 

T := SIMPLIFYSONS(z;,T) 

end for 
else 

test := true 
end if 

if test = true then 
h ■■= Uwew^ 

(R, S) := {N + Pu,N + {h) + P„) {diff-specification associated to (A^, W) in V(P„) } 
t := DIFFTOCANTREE(ii, S) 
hang t from u 

if parent(u) exists and Pparent(«) ^ ^ then 

test := false 
end if 
end if 



Furthermore, GENCANTREE algorithm performs a GCS computation for each list of 
segments whose associated reduced Grobner bases specialize properly, obtaining a tree for 
which the subtrees hanging from the root correspond to the generalized can-specifications of 
the lists configuring a partition of the parameter space. Thus, MCCGS algorithm performs 
the discussion and obtains the Minimal Canonical Comprehensive Grobner System stored 
in the output tree T. 

Example 9. [Singular points of a conic] The general equation of a conic can be reduced 
by a suitable change of variables to the form 



/ 



-|- by^ + 2cxy + dx 



0. 



To study its singular points consider the system of equations 

df df 



S 



5x' dy 



and apply MCCGS algorithm to S using lex(3;, y) and lex(6, c, d) for variables and parame- 
ters respectively. The result is shown in Figure El The interpretation of the output tree is 
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tree T ^ SIMPLIFYSONS(7;, T) 

Input: V a vertex at odd level of tree T where to start the simphfications. 
Input: The tree after simplifications 

Description: 

SIMPLIFYSONS just simphfies the subtree under v on the global T in order to not having 
cancellations nor inclusions between the children of v. Let P be the prime stored in vertex 
V. The simplification is performed as follows: 

Check that there is no Pi child of P such that Pi = Pij . And if any, hang to P all subtrees 
descendant from Pij and drop both Pi and Pij from T. 

Then check whether there is any pair of children ol P, P ^ Pi, P ^ Pj, such that Pi C Pj. 
If so, drop subtree hanging from Pj and also vertex Pj. 

for all V E children(u) do 
if Py = -Pchiid(D) then 

hang from u all subtrees under child(u) 
drop V and child(t;) from T 
end if 
end for 

if there 3 v,w Gchildren(n) such that P^ C then 

drop subtree with root w from T 
end if 



the following. 

There are three different segments: The generic case with Ipp set [1] where the conic has 
no singular points, the segment with Ipp set [y, x] corresponding to a single singular point 
in the conic, and the segment with Ipp set [x] corresponding to a solution with one degree 
of freedom, where the conic is a double line. The conditions over the parameters given by 
the trees are to be interpreted in the following way: 



Ipp 


Basis 


Description 


[1] 


[1] 


\ ((V(6) \ (V(c, b) \ Y{d, c, b))) U V(d)) 




[2cy + d, x] 


(V(6) \ V(c, b)) U (V(d) \ Y{d, b - c^)) 


[x] 


[x + cy] 


V(d,6-c^) 



Figure H] shows the geometrical description of the partition of the parameter space 
provided by the three segments. The generic segment occurs in the whole 3-dimensional 
space except the two planes V(c) and Y{d) plus the line V(c, 6) except the point (0,0,0). 
The one-singular point segment occurs in the two planes V(c) and V(d) except both the 
line Y{c,b) and the parabola Y{d,b — c^). Finally the double line occurs on the parabola 
V(d,6-c2). 



14 



tree t ^ DIFFTOCANTREE((/, J)) 

Input: (I, J) a diff-specification. 

Output: a tree structure containing the Can-Specification of V(/) \ V(J). 



initialize local tree t 

{Pi} :=PRIMEDECOMP(/) 

for all Pi do 

if Pi / VJ + Pi then 

store the Pi as the children of root in t 
{Pij} :=PRIMEDECOMP( J + Pi) 
store the Pij as the children of Pi in t 
end if 
end for 




d 




b 



Figure 4: Geometrical description of the MCCGS for the singular points of a conic system. 
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B{2a, 2b) 










A{0,0) 


jlT ^ s ^ ^ \ \ 

^ ^ ^ r~^^~— — — — _ \I 
^ 1 / 1 ' H§ 






^"^^^^ \ / 1 ^ / \ 
^^\\ ' 1 / \ 

C{2c,2d) 



Figure 5: Nine points circle Theorem. 



5 Applications 

We use now the algorithm to prove part of the 9 points circle Theorem on a triangle. It 
states: For every triangle, the circle through the three middle points of the sides is also 
incident with the height feet. To prove it, and also to obtain supplementary hypotheses if 
needed, consider a triangle with vertices at the points ^(0,0), B{2a,2b) and C{2c,2d) and 
denote P{x,y) the height foot from A (see Figure [5]). The first set of hypotheses are the 
equations of the side BC and the height from A defining the point P{x,y): 

hi: (6 - d)x + (c - a)y + 2ad - 26c = 
/i2 : (c — a)x + {b — d)y = 

Denote r and (xq, yo) the radius and the center of the circle through the three middle points 
(a, 6), (c, d) and {a + c,b + d). Its equation will be {x — xof' + (y — yof' — = 0. So we 
have the three new hypotheses: 

/i3 : (a - xqY + {b- yof - r^ = 
/i4 : (c - xof + id- yof - r2 = 
/i5 : {a + c- xq)^ + {b + d- yo)^ -r'^ = 

The thesis of the theorem is that the circle is incident with the point P{x, y), thus that the 
polynomial 

f = {x- xof + {y- yof - 
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[r2, yO, xO, y, x] 




Figure 6: Generic case for HT in the nine points circle Theorem. 



is zero as a consequence of the hypotheses. The first to do is searching for the solutions of 
the system HT = {hi, /i2 7 ^3; ^4; ^Si /)• Thus we call 

mccgs(ifT, grevlex(x, y, xq, yo, ^2), lex(a, b, c, d)), 

where we set r2 = r^. We obtain a canonical tree with nine cases. But only two cases 
are really interesting. The first one is the generic case (see Figure [6]) for which the Ipp are 
[r2, Ho, xq, y, x] showing that for parameter values not in Y{ad — be) U V((a — c)^ + (6 — c?)^) 
there exists a unique solution. For the real case it is sufficient to consider ad — bc^ 0, as the 
real part of the second variety is inside the first one. The second interesting case is the case 
with basis [1] where no solution exists. The corresponding tree shows that it covers both 
varieties N{ad — be) U V((a — c)^ + (6 — d)^) except for very special cases corresponding to 
degenerate triangles. Thus we have proved that the theorem is true whenever ad — cb ^ 0. 
We can also go further and ask if the thesis is a real consequence of the hypotheses, i.e. if / 
belongs to the radical of the hypotheses ideal H = {hi, h2,hz,hi, /15) whenever ad — bc^ 0. 
To test this we must have 

HTi = {hi,h2, /13, ^4, h, 1 - wf) = (1) 

i.e. the Grobner basis of HTi is [1]. We call now 

mccgs ([/ii, /i2, /13, /14, 1 - wf],grevlex{w,x,y,xo,yo,r), 
lex(a, b, c, d),notnull = {ad — be}), 

and the result is a unique case with basis [1]. Thus effectively / belongs always to the ideal 
of the hypothesis whenever ad — be ^ 0. 

6 Comparison of algorithms 

The CGS of a parametric ideal / can have very different properties as commented in sec- 
tion [H For example 

i) the subsets Si of the parameter space K"^ in which the CGS are divided can be very 
different, they can contain different number of segments, they can overlap, and so on; 

ii) a CGS can contain incompatible segments; 
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iii) the basis Bi can be reduced or not; 



iv) even when a given algorithm does not theoretically ensure some property it can how- 
ever hold it experimentally in most examples. 

So it is quite difficult to make automatic comparisons of the outputs. 

There are three available known implemented methods for obtaining a GCS: 



i) Weispfenning CGB implemented by ( Dolzman-Seidl-Sturm 20061 ) in Reduce. 



ii) Suzuki-Sato SACGB implemented in Risa/AsiiH, 

iii) Montes MCCGS implemented in Maple 8 by M. Manubens in the DPGB library 7.0. 

Even though we use some criteria to evaluate them: correctness of the results, existence 
of incompatibilities, existence of overlaps, number of segments, whether the Si form a 
partition, whether or not specializations preserve the Ipp's of the bases, reduction of the 
bases, theoretical canonicity ensured, theoretical minimality ensured, execution time. 

Although it is not possible to evaluate Weispfenning's CCGB algorithm in practice 
because it has not been implemented, we can analyze its theoretical features. The canonicity 
of CCGB comes from the use of primary decompositions over the conditions, but the method 
is not dichotomic and so the segments are not disjoint. As its objective is to obtain a 
canonical CGB, the bases of the corresponding CGS are faithful and therefore not reduced, 
so specializations do not preserve their 1pp. Furthermore as the segments are not disjoint, 
minimality does not hold. 

For the comparisons with the implemented methods, we have used a Pentium(R) D 
CPU 3.00 GHz, 1.00 GB RAM for the computations and tested different examples using 
the above implementations. 

We have not been able to obtain CGB Reduce in time for these comparisons, so we could 
only test some very simple executions with a demo version. To what we have experimentally 
observed, it gives a partition of the parameter space containing quite more segments than 
MCCGS. The bases are faithful, which is interesting to compute a CGB, but do not give 
direct information on the type of solutions, as these bases are not reduced. It seems to be 
very efficient but the provided results are difficult to be interpreted. In the future we will 
make a more precise analysis. 

SACGB is a very simple and interesting alg orithm based on K alkbrenner's theorem 



for stabilization of polynomial ideals over rings (iKalkbrenner 19971 ) under specialization. 
The published algorithm provides a highly complex CGS, containing even incompatible 
segments, but the Risa/Asir implementation makes an initial reduction and gives a better 
output. We implemented an extra routine to further reduce the output by transforming 
specifications into red-specifications characterized by a pair (N,W), where is the null- 
condition ideal and is a set of irreducible polynomials. 

Among the tests we have done we explain four interesting ones. 

Example 10. First we consider a very simple but illustrative example: the discussion of 
the singular points of a conic already studied in example M 



^There exist also a preliminary Maple version but it is not yet fully developed. 
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Using the Risa/Asir implementation of SACGB together with the additional simplifica- 
tions we obtain the following description of the CGS: 



Ipp's 


Basis 


Description 


[1] 


[1] 


\ W{bcd) 


[1] 


[1] 


V(c) \ Y{d) 


[x] 


[x + cy] 


Y{d,b-c^) 


[y,x] 


[y, x] 


V(d)\V((6-c^)c) 




[y, 2x + d] 


V(d,c)\V(6) 


[y,x] 


[2cy + d, x] 


V(6) \ Y{cd) 



There are two segments with basis [1], i.e. when the conic has no singular points. The 
first one corresponds to the whole space except the three planes V(6), V(c) and Y{d). 
The second one corresponds to the plane V(c) except the line V(c, d). They have empty 
intersection and its union describes the unique generic segment in MCCGS, namely the 
whole space except the two planes V(6) and Y{d) plus the line V(6, c) except the origin 
(0,0,0). 

The segment with Ipp set [x] (i.e. the conic is a double line of singular points) coincides 
with the one in MCCGS. 

Finally, there are three segments with Ipp set [x, y], i.e. the conic has one single singular 
point. The first one corresponds to the plane Y{d) minus the line V(c, d) and the parabola 
V(d, b — c^). The second one corresponds to the line V(c, d) minus the origin (0, 0, 0). The 
third one corresponds to the plane V(6) minus the lines Y(b,c) and Y{b,d). These three 
sets have no common intersection and their union describes the plane Y(b) minus the line 
Y(b, c) plus the plane Y{d) minus the parabola Y{d, b — c^), which is the unique segment in 
MCCGS. Also the basis given by MCCGS for this segment specializes to the bases of the 
three segments provided by SACGB. 

Using Reduce implementation of Weispfenning's CGB, we obtained the following CGS: 



Segment 


Basis 


Description 


1 




b'^cd - bcM ^ 


2 


[x^ + 2cxy + dx + by'^, 2x + 2cy + d, 






cx + by, (26 - 2c'^)y - cd\ 


6 - cV 0, c 7^ 0, 6d = 


3 


[2cdy + d'] 


6/0,d/0,c = 


4 


[x^ + 2cxy + dx + by'^, 2x + 2cy + d,cx + by] 


6^0,c = 0,(i = 


5 


[{2b - 2c^)y - cd] 


c 7^ 0, d 7^ 0, 6 - = 


6 


[x^ + 2cxy + dx + by'^ , 2x + 2cy + d,cx + by] 


c 7^ 0, d / 0, 6 - = 


7 


[icxy + iby'^ - 2cdy - d^] 


(i/0,6 = 0,c = 


8 


[x^ + 2cxy + dx + by'^, 2x + 2cy + d] 


6 = 0,c = 0,d = 



As it can be seen, the description of the segments is not very friendly. In order to 
interpret these CGS as a partition we have manually built the following binary table in 



19 



[x'^3-a, y^A-b, x+y-t] 



[tM2, y, X] 




[f^l 1 , y*t, y'^2, X] [MO, fl^2, y'^2, x] [^6, t'^4*y, y'^2*t'^2, y'^3, x] 



[-4096*b'^3+729*a'^4] [64*b'^3+729*a'^4] 



[b,a] 



[64*b'^3+729*aM] [-4096*b"3+729*aM] [b, a] [b, a] [1] 

Figure 7: Canonical tree for [x^ — a, y'^ — b, x + y — t] wrt lex(x, y, t). 



which represents "being equal to 0", and 1 "being different from 0". The last column 
matches each CGS segment with one of the three MCCGS segments identified by its 1pp. 



Segment 


b 


c 


d 


b-c' 


MCCGS Ipp 


1 


1 


1 


1 


1 


[1] 


2 





1 





1 


[x,y] 







1 


1 


1 






1 


1 





1 




3 


1 





1 


1 


[1] 


4 


1 








1 


[x,y] 


5 


1 


1 


1 





[1] 


6 


1 


1 








[x] 


7 








1 





[1] 


8 














[x] 



The CPU times are 1.46 sec for MCCGS, 0.18 sec for SACGS and 0.05 sec for CGB. 

We see that MCCGS outputs a simpler discussion, not only theoretically but also ex- 
perimentally as all the segments corresponding to the same set of solutions are summarized 
in a single segment, while SACGS and CGB do not. Nevertheless, SACGS and CGB are 
both correct and faster than MCCGS, and although they do not ensure that the Si form a 
partition of the parameter space, in this example they do. 



Example 11. We consider now an example proposed in (jSuzuki-Sato 2006l l for which they 



give the following comprehensive Grobner basis wrt lex{t, x, y) 

S := [x^ — a, y'^ — b, x + y — t] 

and ask for the CGS of (S) wrt lex(j;, y, t). 

MCCGS provides in 632 sec. the canonical tree shown in Figure [7] with only 4 segments 
which takes 16 lines of a Maple worksheet. 
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On the other hand, the Risa/Asir SACGB with the respective simphfications produces 
in 1.62 sec. the following CGS which takes 22 hnes of a Maple worksheet: (for space 
restrictions we do not print the bases) 



Ipp's 


Description 




\ V(a6(729a4 + 646^)(729a4 - 40966^ )(16767a4 + 56326^)) 




V(16767a^ + 56326^) \ V(a6) 




V(a) \ V(6) 




V(6) \ V(a) 




V(729a4 -40966^) \V(a&) 




V(729a4 -646^) \V(a6) 




V(6, a) 



The first segment is described by the whole space minus the three curves V(729a^ + 
6463), ¥(7290^^ - 40966^), V(16767a^ + 563263) and the lines V(a) and V(6). The second 
one is described by the curve V(16767a^ + 56326^) except the origin. The third one is the 
line V(a) minus the origin and the forth segment is described by the line V(6) minus the 
origin. These four segments have empty intersection and are associated to bases with Ipp 
set \t^'^,y,x\. Their union corresponds to the unique generic segment in MCCGS, namely 
the whole space minus the two curves ¥(7290^^ + 6463) ^{7290^ - 40966^). 

The segment in SACGB with Ipp set [t^^^ ,ty,y'^ ,x\ is described by the curve V(729a^ — 
409663) except the origin, which corresponds exactly to the segment associated to the same 
Ipp set in MCCGS. 

The segment with Ipp set [t^'^, t^y, y^, and described by the curve Y{72S)a^ — 6463) 
minus the origin also coincides with the one in MCCGS associated to this Ipp set. 

And finally, the segment having basis with Ipp set [t^, t^y, t^y'^, y^, x] is described on the 
origin V(6, a), which agrees with the segment associated to the same Ipp set in MCCGS. 

All seven segments have no common intersection and thus they form a partition of the 

space, even though SACGB does not ensure it. 

Example 12. We also have tried to test SACGB with the systems of the nine points circle 
theorem explained in section [5] above. SACGB after 3 hours of computation went out of 
memory and had not yet reached an end, while MCCGS takes only 11.45 sec. for testing 
the compatibility of the hypotheses and 2.21 sec. for discussing the theorem thesis. 

Example 13. The last test is the system of the Romin robot ( Gonzalez- Recio 199^ ): 



R= [a + dsi, b - dci, /2C2 + hcs - d, I2S2 + hs^ - c,Si + c^- 1, 53 + C2 - 1, Sg + Cg - 1] 

wrt lex(c3, S3, C2, S2, ci, si) and lex(Z2, ^s, a, 6, c, d). MCCGS takes 43.23 sec in discussing the 
system and provides 9 segments. SACGB also went out of memory. 



Conclusions 

The interest of MCCGS relies, essentially, in the simplicity of the output for applications, 
and in the canonical character of it, conceding an easier interpretation of the results. We 
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have also observed that the obtention of the MCCGS from the BUILDTREE CGS only 
increases the computation time in about 20-30%. 

The existence of the MCCGS depends on the Conjecture formulated in (jMontes 2007I ). 
The use of the algorithm will provide evidence of it or a counterexample. In almost all the 
high number of tests that we have done the algorithm has always obtained a unique segment 
for each different Ipp set, confirming the conjecture. The only ideal for which the algorithm 
obt ains two differ ent segments with the same Ipp is = {u{ux + l),{ux + l)x) proposed 
by (jWibmer 20061 ). and there both segments are clearly intrinsically different and cannot 
be merged nor summarized into a single one. Thus this example also provides evidence of 
the Conjecture. To give a counterexample proving the falsehood of the Conjecture, we must 
find an ideal for which the algorithm MCCGS obtains two or more segments with the same 
Ipp which could be merged or summarized in a different way. 

Although we have only made some very simple tests with CGB, we have observed 
that it seems faster than SACGB and MCCGS in those specific problems. It stands out 
for computing a CGS with faithful bases which are not always useful for applications. 
Experimentally, it seems to obtain a partition of the parameter space, even if there is no 
theoretic evidence. Nevertheless, the number of segments is much higher than MCCGS and 
are difficult to understand. 

SACGB stands out for being in general very reliable to compute a CGS. Its efficiency 
depends on the type of system to be dealt with. It seems to behave faster than MCCGS in 
problems for which a low number of cases is expected. Furthermore, we must remind that 
the output of SACGB is very complex and also needs extra simplifications to be interpreted. 

One can also adapt the MCCGS algorithm to the CGS obtained by other algorithms 
instead of BUILDTREE. To do this one needs to transform the output of the involved 
algorithm into a disjoint reduced CGS, and then apply step ii) and iii), i.e. SELECTCASES 
and MCCGS. 

MCCGS takes, generally, more CPU time for simple problems. Nevertheless the sim- 
plifications inside MCCGS often allow to discuss systems of higher complexity, as seen in 
examples [l2] and [l3] above. 

Finally, we have seen that MCCGS algorithm stands out for having the best features to 
be used for automatic theorem proving and discovering as well as for other applications. 
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